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Abstract. We look at the high- lying eigenstates (from the 10,001st to 
the 13,000th) in the Robnik billiard (defined as a quadratic conformal map of 
the unit disk) with the shape parameter A = 0.15. All the 3,000 eigenstates 
have been numerically calculated and examined in the configuration space 
and in the phase space which - in comparison with the classical phase space 
- enabled a clear cut classification of energy levels into regular and irregular. 
This is the first successful separation of energy levels based on purely dynam- 
ical rather than special geometrical symmetry properties. We calculate the 
fractional measure of regular levels as p\ = 0.365 ± 0.01 which is in remark- 
able agreement with the classical estimate p\ = 0.360 ± 0.001. This finding 
confirms the Percival's (1973) classification scheme, the assumption in Berry- 
Robnik (1984) theory and the rigorous result by Lazutkin (1981,1991). The 
regular levels obey the Poissonian statistics quite well whereas the irregu- 
lar sequence exhibits the fractional power law level repulsion and globally 
Brody-like statistics with (3 = 0.286 ± 0.001. This is due to the strong local- 
ization of irregular eigenstates in the classically chaotic regions. Therefore 
in the entire spectrum we see that the Berry-Robnik regime is not yet fully 
established so that the level spacing distribution is correctly captured by the 
Berry- Robnik-Brody distribution (Prosen and Robnik 1994). 
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1 Introduction 



In the early days of quantum chaos Percival (1973) has suggested, using the 
semiclassical picture and the correspondence principle, that the eigenstates of 
a classically nonintegrable and (partially) chaotic Hamiltonian system should 
be classified as regular or irregular, depending on whether they are associated 
with a classically regular or chaotic regions in the phase space. In fact, he 
referred mainly to the sensitivity of the energy levels with respect to some 
family parameter (of a one-parameter family of Hamiltonians) , by saying 
that in irregular levels the second differences (now known and studied as 
curvature) are typically much larger than in regular levels. This picture has 
been recently confirmed in the studies of curvature distribution of quantum 
spectra of classically nonintegrable systems (Gaspard et al 1990 , Takami 
and Hasegawa 1992) where it has been finally demonstrated that here again 
we encounter universality classes (Zakrzewski and Delande 1993, von Oppen 
1994a,b). 

On the other hand the irregular levels have been supposed to be associated 
with the classically chaotic regions also geometrically in the sense of the Prin- 
ciple of Uniform Semiclassical Condensation (PUSC, see e.g. Li and Robnik 
1994a, Berry 1977a,b), which states that the quantal phase space distribution 
function like Wigner or Husimi or anything in between in the semiclassical 
limit h — > uniformly condenses on the associated classical invariant object, 
which in this case is the underlying chaotic component. According to the 
same principle the regular levels are associated with the quantized invariant 
tori. Thus in the semiclassical limit the entire spectrum and the associated 
eigenstates are decomposed in the regular sequence and generally many irreg- 
ular sequences. More quantatively, it is intuitively obvious to assume that the 
fractional level density of each level sequence is equal to the fractional phase 
space volume of the underlying invariant classical component, and this is one 
of the main assumptions of the Berry-Robnik (1984) theory. It is also a rigor- 
ous result by Lazutkin (1981, 1991) that the torus quantization in the plane 
classical convex billiards works only for states where the classical invariant 
tori exist and the relative measure of such regular states, according to him, 
is equal to the fractional phase space volume of the classically regular regions. 

So far an inspection and separation of spectra in this sense has been per- 
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formed only in one special case (Bohigas, Tomsovic and Ullmo 1990a-b, 1993), 
where the regular levels are characterized by being almost degenerate pairs 
due to some discrete geometrical symmetry. It is the goal of our present 
paper to perform a separation of regular and irregular eigenstates on purely 
dynamical grounds in a generic system of mixed type classical dynamics de- 
scribed by a KAM-type scenario, which is to the best of our knowledge the 
first such analysis. The system that we shall analyze numerically is the Rob- 
nik billiard (1983,1984) (defined as a quadratic conformal map of the unit 
disk) with the shape parameter A = 0.15, and this work is a natural contin- 
uation of our previous work (Li and Robnik 1994b) 

We have calculated the 3,000 eigenstates between the 10,001st and the 13,000th 
of even parity in configuration space and in the phase space and performed 
the classification of states in comparison with the classical dynamics (Poincare 
surface of section plots). After this separation of regular and irregular states 
we have done rather complete spectral statistical analysis which we present 
in section 3. As we shall see the ingredients of the Berry-Robnik picture are 
almost completely confirmed, except for the localization phenomena in the 
chaotic states which are responsible for the deviation from GOE statistics 
to which ideally the statistics must converge in the strict semiclassical limit 
h — > 0. These phenomena have been recently demonstrated, discussed and 
partially (qualitatively) explained by Prosen and Robnik (1994b). 

2 The preliminaries: The calculational tech- 
nique and the method of analysis 

The theoretical questions that we addressed in the introduction are of course 
completely general but in order to illustrate them we have to confine our- 
selves to some specific and possibly generic system. For such a purpose it is 
ideal to study billiard systems. We have chosen the Robnik billiard (defined 
as a quadratic conformal map of the unit disk w = z + \z 2 ) with the shape 
parameter A = 0.15, which is a convex plane billiard. This system and the 
conformal mapping diagonalization technique has been introduced by Robnik 
(1983, 1984) and further studied by Berry and Robnik (1986), Robnik and 
Berry (1986), Hayli et al (1987), Frisk (1990), Bruus and Stone (1994), Stone 
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and Bruus (1993,1994), Markarian (1993), Prosen and Robnik (1993,1994b), 
Li and Robnik (1994a,b) and Backer and Steiner (1994). This one parameter 
system (called Robnik billiard or Robnik model by most workers in the field) 
is a nice generic system which at A = is the integrable circle billiard, for 
0<A<l/4itisa typical KAM system (the KAM theorem applies because 
the boundary is analytic and convex), for A > 1/4 it becomes non-convex and 
eventually for some sufficiently large A < 1/2 becomes ergodic, mixing and 
K-system. Markarian (1993) has recently rigorously proved this property for 
A = 1/2 (cardiod billiard). Li and Robnik (1994c) have numerical evidence 
that ergodicity (and mixing and K property) exists already at A > 0.2775, 
but could set in even earlier. More details about the system and the numer- 
ical technique are given in our recent papers (Li and Robnik 1994a,b). 

In calculating the eigenfunctions in configuration space we have used the 
Heller's method (1991) of plane wave decomposition by applying the sin- 
gular value decomposition (Press et al 1986). All the technical details and 
tricks have been described in (Li and Robnik 1994a,b). The problem with 
missing eigenstates has been overcome by using the exact (double precision: 
16 digits) eigenenergies for all 3,000 consecutive states, from the 10,001st 
to the 13,000th state. The exact energy levels have been obtained by the 
conformal mapping diagonalization technique described in detail in (Prosen 
and Robnik 1993a). 

We have calculated the eigenfunctions in configuration space, their Wigner 
functions in phase space, also their projections onto the surface of section, 
and finally also their Husimi type smoothed objects on the surface of section. 
All the definitions, formulae and the method of presentation are exactly as 
introduced in (Prosen and Robnik 1993b), and later employed in (Li and 
Robnik 1994b). Therefore and because we shall not show any plots of eigen- 
functions in this paper we refer the reader to those previous works. The 
quantal phase space plots (smoothed projections of Wigner functions onto 
SOS) have been compared with the classical SOS plots for all the 3,000 states. 
(Of course we have also plotted the configurational eigenfunctions but this is 
not so important for the classification of states.) One smaller but represen- 
tative part of this output (200 eigenstates) shall be published separately (Li 
and Robnik 1994d) and is also available from the authors upon a request. 
Some representative eigenstates are shown in (Li and Robnik 1994b). 
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In spite of the vast amount of our material (altogether 2x3,000 plots) there 
was no difficulty in classifying the eigenstates. A state is declared regular if it 
is localized on invariant torus belonging to the classical regular region in the 
classical SOS. A state is declared irregular if it is localized (either uniformly 
or non- uniformly) on a classically chaotic region in the SOS. According to 
the PUSC we expect that in the strict semiclassical limit H — > 0, or equiva- 
lently when the energy goes to infinity, the irregular states become uniformly 
extended on the classical chaotic component. However, before this limit is 
reached we observe strong localization of irregular states due to the slow clas- 
sical diffusion inside the underlying classical chaotic component. In fact the 
vast majority of our irregular states are localized chaotic states and only a few 
irregular states are extended chaotic. The consequences of such localization 
for the spectral statistics have been recently discussed and demonstrated by 
Prosen and Robnik (1994b): For the associated irregular level sequence we 
find the phenomenon of the fractional power law level repulsion and globally 
a Brody-like behaviour which then approaches GOE deeper in the semiclas- 
sical limit of smaller and smaller effective h. This is confirmed again in our 
present paper in section 3. 

In performing the classification of states we had 109 cases which at first sight 
would be classified as mixed states in the sense that they "live" both par- 
tially in classically regular regions and partially in irregular regions. However 
a closer examination of finer details of the quantal phase space functions (by 
decreasing the smoothing length) led us to the conclusion that they should 
be classified as regular states. Other a posteriori reasons for such decision 
will be described in the next section. 

The central result of this paper is the counting of regular and irregular states. 
Among 3,000 states we found 987 regular, 109 mixed and 1904 irregular 
states. The percentages are 32.9%, 3.6% and 63.5%, respectively. This is the 
result of the preliminary classification. With the decision mentioned before 
and further explained in section 3 we absorb the mixed states in the set of 
regular states (so now there are 987+109 =1096 of them) which results then 
in pi = 0.365 ± 0.01 for the regular states and p 2 = 1 — p\ = 0.635 ± 0.01 
for irregular states. This is in excellent agreement with the classical rela- 
tive phase space volume of the regular regions as calculated and reported 
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by Prosen and Robnik (1993a) where p\ = 0.360 ± 0.001. We believe that 
this is quite striking confirmation not only of Percival's scheme (1973), but 
also of the main assumption in Berry- Robnik (1984) theory of energy level 
spacings and more specifically of Lazutkin's rigorous results (1981,1991) on 
convex plane billiards. 

In order to get some idea about the convergence of these numerical fig- 
ures with the size of the sample we should note that examination of ad- 
ditional 1,000 consecutive states, namely from the 13,001st to the 14,000th 
state, and the classification by the same criteria as above gives the result 
pi = 0.360 ± 0.01 which is then in brilliant agreement with the classical 
value. However, due to the statistical fluctuations all the statistical spectral 
measures do not necessarily improve monotonically as we shall comment in 
the next section. 

We believe that this is the first detailed dynamical analysis of a large sample 
of eigenstates enabling the separation of regular and irregular states based 
on the classification procedure explained above. This successful separation of 
regular and irregular states makes it possible to analyze the statistical prop- 
erties of regular and irregular sequences separately with the goal to confirm 
the aspects of the Berry-Robnik approach (1984) which is the subject of the 
next section. 

3 Statistics of energy level sequences 

The subject of the Berry-Robnik (1984) theory is the statistics of energy 
spectra of quantal Hamiltonians whose classical counterparts have mixed 
classical dynamics in the sense that classical regular regions and classical 
irregular regions coexist in the phase space. It is based on our knowledge 
of spectral fluctuations and their statistics in the context of quantum chaos 
(Berry 1983, Giannoni et al 1991, Haake 1991, Gutzwiller 1990, Eckhardt 
1988, Bohigas and Giannoni 1991, and Robnik 1994). We know that there 
are three universality classes of spectral fluctuations: Poisson statistics in 
the classically integrable cases; in case of classical ergodicity we find the 
GOE/GUE statistics of random matrix theories depending on whether there 
is one/none antiunitary symmetry (we ignore spin). The interesting and dif- 
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ficult case of mixed type classical dynamics of KAM-like (generic) systems 
has been studied for the first time by Robnik (1984) numerically, where a 
continuous transition from Poisson to GOE statistics in a billiard system 
(Robnik 1983) has been found, and this work has been substantially revised 
in (Prosen and Robnik 1993a). Further theoretical progress was published 
by Berry and Robnik (1984) where the following semiclassical theory of the 
level spacings has been presented. The eigenstates (their Wigner functions 
in phase space) are supposed to condense uniformly on the underlying clas- 
sical invariant regions such that each of them — in the semiclassical limit - 
supports a level sequence which for itself has a Poisson or GOE statistics if 
the region is regular or irregular, respectively. All the regular regions can be 
thought of as supporting a single Poisson sequence because the Poisson statis- 
tics is preserved upon a statistically independent superposition. The mean 
level spacing of such a sequence is determined by the fractional phase space 
volume of the regular regions. On the other hand each chaotic (GOE) level 
sequence has a mean level spacing governed by the corresponding fractional 
phase space volume. The entire spectrum is then assumed to be a statistically 
independent superposition of all subsequences. The statistical independence 
in the semiclassical limit is justified by the principle of uniform semiclassi- 
cal condensation of eigenstates (in the phase space) and by the lack of their 
mutual overlap, in consistency with Percival's (1973) conjecture. Thus the 
problem of the statistics of the entire spectrum is now mathematically pre- 
cisely formulated (this forms the essence of the Berry-Robnik approach) and 
its solution as far as the level spacings are concerned can be expressed in 
the following way: The statistical independence of superposition implies fac- 
torization of the gap distribution functions (Mehta 1991, Haake 1991): The 
probability that there is no level within a gap clearly factorizes upon a statis- 
tically independent superposition. The connection between the level spacing 
distribution P(S) and the gap distribution E(S) is as follows 



and conversely 



Leaving aside the general case of infinitely many chaotic components which 
does not include anything surprisingly new let us restrict to the case of one 




(2) 
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regular component with mean level density p\ (= fractional phase space 
volume) and one chaotic component with the mean level density p 2 where 
Pi + P2 — 1- This is going to be already an excellent approximation because 
in a generic system of a mixed type there is usually only one large and 
dominating chaotic region. Following (Mehta 1991, Haake 1991 and Berry 
and Robnik 1984) we have 

E(S) = E Poisson ( Pl S)E GOE (p 2 S) (3) 

where the Poissonian gap distribution -Ep isson is 

£poisson(S) = exp(-S) (4) 

whereas for the Eqqe there is no simple closed formula (for the infinitely 
dimensional GOE case) and must be worked out by using practical approx- 
imations for Pgoe and/or Eqoe which e.g. can be found in (Haake 1991, 
pp72-74). However 2-dim GOE case (the so called Wigner surmise) can be 
worked out explicitly as given in (Berry and Robnik 1984, formula (28)) 
which usually is a good starting approximation. 

As for the delta statistics A(L) the similar procedure based on the as- 
sumption of statistical independence leads to the simple (additive) formula 
(Seligman and Verbaarschot 1985) 

A(L) = Ap oisson ( Pl L) + A GOE {p 2 L) (5) 

where A Poisson (L) = L/15 whilst for A GO e there are good approximations 
given in (Bohigas 1991). 

The main objective of this paper is to verify the elements and aspects of 
the Berry-Robnik approach. However, before such a regime is formed in a 
KAM system in the ultimate far semiclassical limit we typically observe a 
quasi- universal behaviour in the spectral statistics which is characterized by 
the fractional power law level repulsion and globally the adequacy of Brody 
(1973,1981) distribution and of similar distributions such as Izrailev (1989). 
A thorough numerical study of this phenomenon has been published recently 
by Prosen and Robnik (1993a, 1994a,b), and will be discussed later on. In 
cases where the chaotic component becomes large (p 2 ~ 1) and the quantal 
chaotic states strongly localized the phenomenological Berry-Robnik-Brody 
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distribution proposed in (Prosen and Robnik 1994b) can be very efficient in 
capturing the global features of the experimental and numerical data (Lopac 
et al 1992, 1994). 

As explained in section 2 we have separated 1096 regular levels from the 
remaining 1904 irregular levels. Now we perform the statistical analysis of 
each subsequence in the next two subsections 3.1-3.2 and then also of the 
entire spectrum in subsection 3.3. 

3.1 Regular levels 

Of course, before performing the separation procedure we have carefully un- 
folded the entire spectrum of 3,000 levels by using the Weyl formula (with 
perimeter and curvature corrections) for even parity states given in formula 
(4) of our paper (Li and Robnik 1994a). After extracting 1096 regular lev- 
els and calculating the (nearest neighbour) level spacings we have first to 
normalize the empirical level spacing distribution P(S) to the unit first mo- 
ment, i.e. such that the mean level spacing is unity. By doing this we get 
again the relative level density of the regular component with the value of 
Pi = 0.365. Having performed this normalization we calculate the cumula- 
tive level spacing distribution W(S) = f$ P(x)dx and compare it with the 
Poisson statistics and also try to fit it with the best fit Brody distribution 

W B (S, P) = l- exp(-6^ +1 ), b = {r((/3 + 2)/(/3 + (6) 

and best fit Berry- Robnik distribution (|3|). The result is shown in figure 
1, where we also show the so-called U-function introduced in (Prosen and 
Robnik 1993a), defined as 

U(W) = — arccos Vl-W, (7) 

7T 

in reference with the best fit Berry-Robnik distribution. As we see the agree- 
ment of our numerical data with the Poissonian distribution is reasonably 
good in the sense that Poisson is within the Ha statistical fluctuations of 
the data. The best fit Brody distribution and the best fit Berry-Robnik dis- 
tribution deviate slightly from Poisson: For Brody we get (3 = 0.041 ±0.002, 
which ideally should be zero, and for Berry-Robnik we get p\ = 0.704, which 
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ideally should be unity. In figure 2 we show the delta statistics for the 
same sequence of 1096 regular levels and for sufficiently small L < L max (for 
L larger than L max the saturation effects set in (Berry 1985): in our case 
Lmax ~ 10) we see excellent agreement with Poisson. The best fit Seligman- 
Verbaarschot delta statistics @ gives p\ = 0.9998, which is surprisingly close 
to the ideal value unity. Thus we think that our statistical analysis of the 
regular sequence of levels clearly supports the theoretical assumptions in the 
Berry- Robnik (1984) approach. 

3.2 Irregular levels 

In analogy to the previous procedure we have extracted 1904 irregular levels 
and normalized the level spacing distribution to the unit first moment which 
yields p2 = 0.634 ± 0.002 in consistency with pi = 0.365 ± 0.002, such that 
we have unit level density of the entire spectrum, namely p± + p2 = 1- The 
cumulative level spacing distribution for the irregular sequence is plotted in 
figure 3, where we see that it deviates from the ideal GOE distribution sub- 
stantially and significantly. This is a consequence of the strong localization 
of irregular eigenstates on classically chaotic components demonstrated for 
even higher energies (ten times larger) in (Li and Robnik 1994b) for the same 
billiard with A = 0.15, whereas the implications for level statistics have been 
demonstrated and discussed by Prosen and Robnik (1994b). The best fit 
Berry-Robnik distribution yields p\ = 0.374 which however ideally should 
be zero. On the other hand the more interesting best fit Brody distribution 
captures the fractional power law level repulsion with (3 = 0.286 ±0.001, and 
as we see both in the W(S) as well as in the U-function plot it is globally 
much better fit than Berry-Robnik which is qualitatively well understood in 
(Prosen and Robnik 1994b). If we were able to go to much higher energies, 
say 100 times higher, then the statistics of irregular levels is predicted to 
approach GOE statistics, i.e. the Brody parameter (3 goes to one. 

As for the delta statistics the analysis of the data by the best fit Seligman- 
Verbaarschot formula fl5|) gives p\ = 0.373, which is equal to and consistent 
with the Berry-Robnik fit for level spacings, but of course is at variance with 
the ideal value p\ = 0. This is certainly the consequence of the localization 
effects. Having understood the reasons for the deviation of spectral statistics 
from the limiting (as Ti —>■ or equivalently as the energy goes to infinity) 
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GOE statistics we may conclude that the relevant aspect of the Berry-Robnik 
(1984) theory is reconfirmed supporting our claim that this theory is asymp- 
totically exact theory. 

3.3 The entire spectrum 

Ideally in the strict semiclassical limit of sufficiently high energies we would 
expect the statistics of the entire spectrum of 3,000 levels (regular plus irreg- 
ular) to be described by the Berry-Robnik distribution (§). This prediction 
is based on two facts: (i) GOE applies for irregular levels, and (ii) regu- 
lar and irregular levels are superposed statistically independently. Now, the 
first assumption is not satisfied due to the localization effects explained in 
the preceeding subsection, so that instead of GOE we find in fact Brody 
with (3 = 0.286 ± 0.001, as shown in figure 3. The second assumption of 
statistical independence is accepted as valid and thus the reasoning leading 
to the factorization of the gap distribution is valid except that we find now 
not Berry-Robnik distribution but the so-called Berry- Robnik-Brody (Prosen 
and Robnik 1994b) distribution instead 

E BRB (S, Pl ,P) = E Pmsson (piS)E B (p 2 S,P) = e-^ s E B ((l- Pl ))S,P), (8) 

where Brody statistics has gap distribution which can be expressed in terms 
of incomplete Gamma function Q 

™=«(^.( r o<r)- <■» 

Indeed this phenomenological two-parameter distribution provides an excel- 
lent fit to our spectral data of 3,000 consecutive levels as can be seen in 
figure 5, with the best fit parameter values p\ = 0.309 and j3 = 0.370. They 
differ slightly from the ideal values p\ = 0.365 and (3 = 0.286 ± 0.001, but 
are consistent with them thereby confirming the statistical independence as- 
sumption. 

Nevertheless, just as a consistency check we can try to fit the spectrum 
with the Brody distribution and the Berry-Robnik distribution. The result 
is shown in figure 6 and the best fit parameters are j3 = 0.166 ± 0.007 and 
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pi = 0.504, respectively. The fits are not bad and Brody is certainly better 
than Berry- Robnik as can be seen especially in the U-function plot. The 
conclusion is that as a consequence of the localization effects the entire spec- 
trum is indeed best described by the Berry- Robnik-Brody distribution fitted 
in figure 5. The value of x 2 i n the latter is about three times smaller than in 
Brody and six times smaller than in Berry-Robnik. 

Finally we should mention that the best fit Seligman-Verbaarschot (|5|) delta 
statistics fitted to the spectrum within the interval 0.9 < L < 9.9, below the 
saturation region, gives p\ = 0.472 which is consistent with the Berry-Robnik 
fit with pi = 0.504 as explained above. For the sake of completeness we show 
the delta statistics in figure 7. 

Unlike the counting measure of regular levels the statistical measures of the 
spectral subsequences do not converge so uniformly and so fast. Calculating 
the level spacing distribution and the delta statistics of an enlarged sample 
(namely 1,000 more levels added on top to the block of 3,000 levels) led ac- 
tually to slightly worse results. 

Finally, as promised, we would like to give additional a posteriori reasons 
for classifying the preliminary mixed type states as regular. (1) Due to the 
strong localization any ambiguous mixed states is likely at least to mimic a 
regular state. (2) After absorbing the mixed states into the set of regular 
ones the relative fraction of regular states pi = 0.365 agrees much better 
with the classical estimate pi = 0.360. (3) The statistical measures (P{S) 
and A(L)) also agree then much better with Poisson statistics. (4) The sta- 
tistical measures of irregular sequence agree then notably better with the 
trend towards GOE in the sense that the Brody exponent (3 increases from 
P = 0.276 to = 0.286. 

4 Discussions and conclusions 

The main goal of our present paper is to analyze a large sample of eigenstates 
of an autonomous Hamiltonian system sufficiently high in the semiclassical 
limit and to perform the classification of the states into regular and irregular 
by examining their Husimi-type phase space plots in correspondence with the 
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classical surface of section plots. This plan has materialized in the specific 
choice of the Robnik billiard with the shape parameter A = 0.15, which is a 
KAM-type system (the boundary is convex and analytic) with one large and 
dominating chaotic component with the relative phase space volume (not 
SOS area) p 2 — 0.640. We have plotted 3,000 consecutive states (from the 
10,001st to the 13,000th) in configuration and in phase space and performed 
the said classification with the result that the relative fraction of regular 
states is equal to 0.365 which is in excellent agreement with the classical 
value pi = 0.360 = 1 — p^. Further we have calculated the statistical mea- 
sures (level spacing distribution and the delta statistics) for the regular and 
irregular sequence separately and also of the entire spectrum. The aspects 
and assumptions of the Berry-Robnik (1984) approach have been confirmed 
although we are not yet far enough in the semiclassical limit to see the GOE 
statistics for the irregular levels, whilst the regular levels are seen to obey 
the Poisson statistics very well. We believe that this is the first direct dy- 
namical separation of regular and irregular levels in a generic system and the 
first such statistical analysis. It confirms the Percival's (1973) classification 
scheme, the ingredients in the Berry-Robnik theory (1984) and the specific 
rigorous results on convex billiards by Lazutkin (1981,1991). It supports our 
claim that Berry-Robnik (1984) theory is asymptotically exact theory. One of 
the future projects should be in going further into the semiclassical limit for 
which new methods are needed and one such might be the employment of the 
quantum Poincare mapping (Bogomolny 1992, Schanz and Smilansky 1994, 
Prosen 1994a,b,c,d). This might lead us to a direct demonstration of the 
applicability of the Berry-Robnik theory with high numerical and statistical 
significance which so far was successful only in an abstract time-dependent 
Hamiltonian system, namely the compactified standard map (Prosen and 
Robnik 1994a,b). 
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Figure captions 



Figure 1: The cumulative level spacing distribution (a) and U-function (b) 
of 1096 regular levels. In (a) the dotted curves represent the Poisson and 
GOE statistics. The thick curve is the numerical data, while the thin curve 
shows the best fit Berry-Robnik distribution and the broken curve shows 
the best fit Brody distribution. The best fit parameters are p\ = 0.704, 
P = 0.041 ± 0.002, respectively. In (b) we display the U-function difference 
U(W) - U(W BR ) versus W(S), where U br = U(W BR (S)) is the U-function 
of the best fit Berry-Robnik distribution. The numerical noisy curves (the 
average value with ±la band) in (b) are compared with the best fit Brody 
distribution (broken curve) and the Poisson distribution (dotted curve). 

Figure 2: The delta statistics of 1096 regular levels and the best fit Seligman- 
Verbaarschot formula (||). The dotted curves are the Poisson and GOE, the 
thick curve shows the numerical data, the best fit Seligman-Verbaarschot for- 
mula is represented by the thin curve which overlaps the limiting Poissonian 
at L < 7.5. The vertical dashed lines indicate the region where the least 
square fit has been done. The best fit parameter value is pi = 0.9998 which 
is in excellent agreement with the ideal value 1. 

Figure 3: The same as figure 1 but for the 1904 irregular levels. The best fit 
parameter values are p\ = 0.374 and (3 = 0.286 ±0.001 for the Berry-Robnik 
and Brody distributions, respectively. In the U-function plot (b) we see that 
the Brody fit is globally very good. (We do not show the Poissonian curve 
here.) 

Figure 4: The same as figure 2 but for the 1904 irregular levels. The best 
fit parameter value is p\ = 0.373 which is at variance with the ideal value 0, 
but is completely consistent with the Berry-Robnik fit in figure 3. 

Figure 5: The cumulative level spacing distribution (a) and U-function plot 
(b) of the total 3,000 regular and irregular levels compared with the best 
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fit Berry- Robnik-Brody formula (P). In (a) the dotted curves represent the 
Poisson and GOE statistics. The thick curve is the numerical data, while 
the thin curve shows the best fit Berry- Robnik-Brody distribution. The best 
fit parameters are p\ = 0.309 and (3 = 0.370. In (b) we show the U-function 
difference U(W) - U(W BRB ) versus W(S), where U brb = U(W BRB (S)) is the 
U-function of the Berry- Robnik-Brody. The thick curve represents numerical 
data and the noisy curves indicate the ±lcr band. The quality of the fit is 
really excellent. However, the ideal values of p\ and (3 would be pi = 0.360 
and, according to figure 3, (3 = 0.286. 

Figure 6: The same as figure 3 but for the total 3,000 levels. The best fit 
parameter values are p\ = 0.504, (3 = 0.166 ± 0.007. 

Figure 7: The same as figure 2 but for the total 3,000 levels. The best 
fit parameter value is p\ = 0.472, which is in complete agreement with the 
Berry- Robnik fit in figure 6. 
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